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Abstract 



Excitons confined to flat semiconductor quantum dots with elliptical cross 
section are considered as we study geometrical effects on exciton binding en- 
ergy, electron-hole separation, and the resulting linear optical properties. We 
use numerical matrix diagonalization techniques with appropriately large and 
optimized basis sets in an effective-mass Hamiltonian approach. The linear 
optical susceptibilities of GaAs and InAs dots for several different size ratios 
are discussed and compared to experimental photoluminescence spectra ob- 
tained on GaAs/ALjGai-zAs and InAs/GaAs quantum dots. For quantum 
dots with several nm in size, there is a strong blue shift of the luminescence 
due to geometrical confinement effects. Also, transition peaks are split and 
shifted towards higher energy, in comparison with dots with circular cross 
sections. 
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I. INTRODUCTION 



The development of clever fabrication techniques in semiconductors has brought the re- 
duction of the effective dimension of electronic states from their usual three-dimensional 
character in bulk materials, to "zero-dimensional" states in quantum dots.0 The quantum 
effects of these lower- dimensional systems have attracted much attention in recent years, due 
in part to possible applications which include electronic devices based on parallel and perpen- 
dicular transport, quantum well lasers, and optical devices.! Two-dimensional quantum- well 
or quantum-film structures, which provide confinement in one space dimension, have been 
well investigated, and quantum exciton effects observable even at room temperature have 
been studied.i The confinement of excitons has also been shown to result in very large 
electro-optical shifts of the absorption peaks, producing the so-called quantum-confined 
Stark effect. 10 

In quasi-zero-dimensional quantum dot systems, the additional quantum confinement 
dramatically changes the optical and electronic properties, compared to those in higher- 
dimensional structures, as the whole single-particle spectrum is now discrete. Correspond- 
ingly, the excitonic spectrum is expected to be strongly affected. The properties of excitons 
confined in quantum boxes were first analyzed theoretically by Bryant,! who used variational 
and configuration-interaction representations. Later on, excitons and biexcitons have been 
studiedjii as well as excitons in the presence of a strong magnetic field,! using numerical 
matrix diagonalization schemes. 

On the experimental side, interband optical spectroscopies, such as photoluminescence, 
have been used to study various quantum dot systems — such as those produced in the 
GaAs/Al x Gai_ x As structure, with its bandgap modulation. More recently, fascinating 
studies on so-called "self-assembled" quantum dots, such as InAs and In^Gai-^As clusters 
on GaAs substrates, have also been reported.EUll Most of the theoretical investigations are 
based on the assumption that the shape of quantum dots is a simple sphere or box, having a 
great deal of symmetry, both because it simplifies calculations and because quantities such 
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as the exciton binding energy scale very well with the overall dot size. However, realistic 
dot shapes are probably much less symmetrical, as well as being typically flat and more 
two-dimensional in shape.&~E! 

Here, we consider the effect that less symmetric structures, namely flat quantum dots 
with elliptical cross sections, or "elliptical quantum disks", have on the excitonic optical 
properties. To date, little work has been reported on the properties of nonsymmetric quan- 
tum dots, probably because this system has more complicated solutions. 00 Our studies 
within the effective mass approximation yield some very interesting consequences of the 
elliptical asymmetry: apart from the expected blue shift of the first excitonic transition 
for dots with the same overall area but different axes, we find a rearrangement of the os- 
cillator strength which characterizes individual dot shapes. In particular, since elliptical 
cross-section dots have less symmetry, some of the accidental degeneracies in circular dots 
giving rise to stronger and fewer peaks in the imaginary part of the optical susceptibility 
are split. This gives rise to a more monotonically decreasing peak intensity for higher en- 
ergy features in the susceptibility of noncircular dots. This behavior can in turn be used to 
structurally characterize specific dots from their photoluminescence excitation response. 

The remainder of the paper is organized as follows. We introduce the theoretical method 
in Sec. II. Here, we outline the effective mass Hamiltonian approach and introduce the 
various basis function representations which allow us to use numerical methods to calculate 
the eigenvalues and eigenfunctions of excitons in these quantum dots. A great deal of care 
is needed to assure that the solutions obtained are well behaved and converged with a finite 
computational effort. We discuss in this section how this is accomplished. In Sec. Ill, we 
discuss the main geometrical effects on various exciton characteristics, such as the exciton 
binding energy, electron-hole separation, and the linear optical susceptibility. Solutions for 
excitons in quantum dots with both circular and elliptical cross sections are shown, using 
large enough basis sets, and a set of optimized basis functions, which improve the accuracy 
of the solutions at a modest computational cost. Finally, we summarize our conclusions in 
Sec. IV. The Appendix contains an outline of the derivation of the Coulomb matrix element 
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with these basis functions. The analytical expression presented there greatly simplifies our 
calculations. 



II. THEORETICAL METHOD 

For concreteness, and to simulate recent quantum dot systems,0"0 we assume quantum 
dots with an oblate spheroidal profile where the lateral xy confinement is much weaker (or 
larger size) than that along the z direction. Correspondingly, the electrons and holes are 
assumed confined in an effectively two-dimensional potential with a constant z profile, V z . 
We assume V z to be a hard-wall confinement potential, so that the z component of the energy 
is h 2 ir 2 /2mL 2 z) with L x , L y 3> L z . Further, we approximate the single z wavefunction in the 
problem as a 5 function centered at the origin, so that the problem can be described by a 
separable Hamiltonian in two dimensions. The lateral confinement is modeled via harmonic 
potentials with two different frequencies uj x and u y , which yield the elliptical cross section 



of the dots with axes ratio given by L x /L y = yJco y /co x , for both electrons and holes. The 
smoothly varying potential should mimic well the situation in experiments where the dots 
are effectively embedded in a dielectric matrix.0lll 

The effective-mass parabolic-band Hamiltonian for an electron-hole pair is given by H = 
H e + H h + H e _ h , where the subscripts e and h represent electron and hole, and 



u P i 1 2 2,1 22,i- i -i \ 



with a similar expression for the Hamiltonian of the hole, The Coulomb interaction 

between electron and hole is screened by a background dielectric constant e, so that H e ^h — 
-e 2 /er e _ h . 

We rewrite the Hamiltonian into relative and center of mass coordinates, described by 
r = r e — Th, and R = (m e r e + m^r/J/M. The total and reduced masses are given by 
M = m e + rrih, and \x = m e mh/M, respectively. The total Hamiltonian of this system can 
then be written in the form H = H c m + H reh with the expected expressions 
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H c . m . = — + -Mu 2 x X 2 + -Mu; 2 y Y 2 + V z , (2) 



2M 2 2 
and 

#re« = |- + T^K^ + / \ = g + ^ 2 • 

III I 2 y ev^c + ?r 

The Hamiltonian of the center of mass is obviously a two-dimensional harmonic oscillator 
in the plane, with wavefunction ^n x n y = 4>n x {.X) 4>N Y (X)y an d energy E cm , where 

MX) = (^m) V2 e~^H N (a M X) , (4) 



&m = \/Muj x /h, and 



E c . m . = (N x + l) hu x + (n y + fuv v + ■ (5) 
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Here, iVx and iVy are quantum numbers for the center of mass coordinate, and is a 
Hermite polynomial.^ 

The physics of the problem is determined to a great extent by the ratio between the 



effective Bohr radius, a* B = h e/fie 2 , and the size of the dot, L = JL x L y , where L 



h/fiUi. The strong confinement limit for L < a* B is characterized by a weak electron- 
hole correlation and by the Coulomb term being a small perturbation of the single-particle 
confined-level energy On the other hand, the weak- confinement limit for L > a* B reduces 
asymptotically to the problem of a free two-dimensional exciton for large L, where the 
Coulomb interaction dominates the state of the excitoJ. 

With this in mind, the effects of the Coulomb term Hi in H re i = H + Hi, are treated by 
using the solutions of H as the basis set in the diagonalization of H re [. The unperturbed 
Hamiltonian of the relative coordinate H is also a two-dimensional harmonic oscillator, so 
that the wavefunction of the interacting electron-hole pair is described by a linear combina- 
tion of wavefunctions, ip nx n y = 4>n x {x)(f>ny(y), with the 0's satisfying a similar expression to 
Eq. (4), and correspondingly 

El, = [n x + i) hu x + [n y + I) fuj v + , (6) 
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where n x and n y are quantum numbers for the relative coordinate. (Notice z confinement 
length for this coordinate is 2L Z .) 

With this basis set, the interaction matrix elements of the electron-hole pair can be 
calculated analytically and expressed in terms of hypergeometric functions as outlined in 
the Appendix. This analytical expression greatly simplifies the calculation, as most of the 
computational time is spent on the calculation of the matrix elements rather than on the 
diagonalization of the matrix. The resulting Hamiltonian matrix is real, symmetric, and 
sparse. The energies and eigenf unctions are calculated from the numerical diagonalization 
of the matrix, for a given size of the basis. The diagonalization is repeated with larger basis 
sets until the desired convergence is achieved (see below). 



A. Circular limit 

As an additional test of our numerical procedures we use the resulting radial equation of 
the relative Hamiltonian for the circular dot case, u x = u y = u, which is then given by 

The resulting one- dimensional equation can be directly integrated numerically, as reported 
by Que.0 Furthermore, the problem can also be solved using a harmonic basis set of radial 
states as those described above. For a large enough basis set, this approach yields the same 
results as those obtained by direct numerical integration^ Notice also that the large dot limit 
(u ~ 1/L ~ 0) is easily solved as an expansion in terms of the Laguerre polynomial-based 
solutions of the free exciton problem.! This allows one to obtain very accurate solutions for 
1/L « with little computation. These results, moreover, allow us to test the convergence 
of the irreducible two-dimensional nonsymmetric problem {u x ^ u y ), as we discuss below. 
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B. Optimized basis 



Solutions for these elliptic cylinderlike quantum dots are also carried out using an opti- 
mized basis set, as the solution method discussed above requires a rather large number of 
states for full convergence, especially for large dots. Here, one notices that for large dots it is 
the Coulomb interaction that dominates the exciton states (since confinement becomes less 
important). Correspondingly, we choose a set of optimized frequencies, Q x and Q y , which 
are larger than the original frequencies, u x and u y . The values of the Si's are determined 
variationally, and allow one to consider if-matrix systems that are much smaller than those 
required when one uses the u basis. The physical reason for this is that as the dot size 
increases, the exciton size converges to a 2 ^ (the radius of the free two-dimensional exciton), 
and one needs a large number of uj states to describe the small-scale structure of the exciton. 

Notice that the harmonic f2 basis allows also an easy calculation of the if -matrix elements 
in this case, so that for example, 



(njriy \H \n x n y ) n 

h 



hQ x (n x + i) + hQ y (n y + -) 



1 



h 



(ni - ui)(n x + -)- — (n; - ^)( % + -) 



„'>»» 



n x (n x - 1) 5 nx/:nx -2 + \lin x + l)(n x + 2) 5 KtJlx+2 

{Q 2 y - U 2 y ) \jn y {n y - 1) 5n yl ,n y -2 + V / K + 1 )K + 2 ) 6 n y/ ,n y +2 



$n y ,,n y 



(8) 



Notice this reduces to the obvious diagonal matrix for Q — > uj. This expression allows one 
to evaluate the Hamiltonian matrix rather conveniently, even for this other basis set. 



C. Exciton characteristics 

The wavefunctions of the relative coordinate problem can then be written as = 
J2n x n y a n x n y \n x ,n y ), with either the u or Q basis states, which can be used to study var- 
ious characteristic properties of the exciton system. For example, the mean electron-hole 
separation r s , is given by 
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which gives an idea of the exciton size. 

One can also use the diagonalization results to calculate directly measurable properties, 
such as the linear optical susceptibility of the quantum dot/disk. The linear optical suscep- 
tibility is proportional to the dipole matrix elements between one electron-hole pair j state 
and the vacuum state, (0|P|l)j.i These in turn are proportional to the interband matrix 
element, Pcuj^ which is the matrix element formed between an electron and hole in the 
conduction and valence bands, respectively. The form of the dipole matrix elements for a 
single exciton in the envelope function approximation is given bj 



|(o|P|i>| 2 = b w | 2 |Vi(o)| s 



J J^(X e ,Y e )dX e dY e 



(10) 



Here, the wavefunction for the relative coordinate is given as above, so that 



|^(0)| 2 = (n/hir)y/w x w y 



n x n y 



(11) 



where the a Uxn coefficients are obtained from the diagonalization of the relative-coordinate 
Hamiltonian, and 



$>{X e ,Y e )dX e dY e 



^ 2 hN x \N Y \ [ 7 rM v ^^2^ + ^(iV x /2)! 2 (iV y /2)! 2 l \ (12) 



with Nx = even and Ny = even, for nonzero matrix elements. Finally, the dipole matrix 
elements have the form 



12 V> 



|<0|P|1>| 2 = 4|^| -iv;. - 



N X \N V \ 



2 n x +n y (N x /2)\ 2 (N Y /2y/ 



(2 x y T7' X \ Tlyl) I a n x n y 



(13) 



The linear optical susceptibility can then be calculated from 



X (co) = J2 |(0|P|l)/(^ - E j - ihT)- 1 , (14) 
j 

where T is introduced as a phenomenological level broadening constant. 

Notice also that since experimental systems are typically configured to analyze a large 
collection of nearby dots, one should in principle be concerned by the effect of local fields. 
However, in typical experimental systems so far, where the separation between dots can 
be several microns, it is valid to assume that dots are basically independent. In the case 
of higher-dot densities, however, the dynamical response of the system may be affected by 
the local fields produced by neighboring dots, and one can obtain that response from the 
individual microscopic polarizabilities.0 



III. RESULTS 

As an interesting example of a typical system, we use parameters to describe GaAs 
quantum dots, so that the dielectric constant is e = 13.1, and the carrier masses are m e = 
0.067m , and (the heavy hole effective mass) m hh = 0.37m E We present the numerical 
results for heavy-hole excitons in GaAs quantum dots with elliptical and circular cross 
sections. The solutions can be calculated using a sufficiently large basis set and/or an 
optimized basis set, as described above. Results obtained from these methods and different 
basis sets are shown in Fig. 1. The exciton binding energy and normalized electron-hole 
separation are shown as a function of quantum dot size ranging from 2 to 100 nm. 

In the insets, results are shown for circular dots, with dot, dashed, and dot-dashed curves 
showing results for basis sets with M = 30, 100, and 500 wavef unctions, respectively. Here, 
states with n and n' from to 29, 99, and 499 are used in Eq. (]?]). (The matrix size is 
obviously M x M, and is diagonalized by a QL decomposition technique.^) The results of 
the one-dimensional radial equation in the weak-confinement limit are shown with the solid 
line for comparison, and represent the exact quantity (both Ei, and r s ) for large L. The 
transition between the strong and weak confinement regimes comes appropriately when the 
size of the quantum dot is near the effective Bohr radius, L « a* B = 12.2 nm. Notice that 
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it is for L pa 15 nm that the M = 100 curve (dashed) departs from the exact result (solid 
line), and that one requires larger M values as L increases to achieve better convergence. 
The M = 500 basis set (dot-dashed curve) yields the convergent solutions with acceptable 
accuracy and execution time for a larger range of L values (< 60 nm). Similar behavior can 
be seen in the electron-hole separation in the ground state [inset in panel (b)]. 

The inset in 1(a) also shows the difference between the uj basis and Q basis results. 
Diamonds show results for the M = 400 basis set with optimized Q frequency, while triangles 
show results for M = 400 with the original u basis set. For both cases, the states with n x , 
n y , n x ', and n y ' from to 19, respectively, are included for the M = 400 basis set (see Eq. 
H). These states are used as the basis set for elliptical quantum dots. The results with 
the optimized basis set are essentially identical to the convergent solutions. The optimized 
frequencies used are also given as a function of dot size in the inset 1(a) (with Eq = hQ) as 
a long-dashed curve. For dot sizes 30 nm and larger, the optimized frequency Q converges 



to 80 meV, corresponding to a dot size Lq = ^Jh/fifl = 4.1 nm. This latter value is close to 
the two-dimensional effective Bohr radius (a^ D = a* B /2\/2 = 4.3 nm), as one would expect. 
Notice that for most of the range of L's shown, one is in the weak-confinement regime, 
where L > a 2 ^ . This is why the optimized basis set gives very reasonable results with 
minimal effort. It is also interesting that the agreement continues also for smaller dot sizes, 
completing the range from 2 to 100 nm. The exciton ground state binding energies and 
the normalized electron-hole separation obtained with the optimized Q basis approach are 
basically exact to the fully converged results. 

The main panels in Fig. 1 (a) and (b) show the geometrical confinement effects of excitons 
with several different size axis ratios in the xy plane. The plots show results versus L = 



^L x L y , the effective size of the dot, for u} x /uj y = 1, 4, and 9 (L y /L x = 1, 2, and 3) with 

diamonds, pluses, and triangles, respectively. The exciton binding energy increases for a 
small L as the axis ratio is increased; meanwhile, the normalized electron-hole separation 
is basically unchanged, except for small L values, where the confinement energy dominates. 
Notice that as the axis ratio increases, the single-particle and exciton states move up in 
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energy but the binding energy increases. This increase in binding energy for the elliptical 
dots is then related to the increase of the Coulomb energy relative to the confinement 
contribution to the total energy. (In fact, since r s /L is basically unchanged with geometry, 
the Coulomb interaction energy is nearly constant in all these cases.) 

To better explore the geometrical confinement effects on excitons, we show the linear 
optical susceptibility of the GaAs quantum disks with elliptical cross sections and lateral 
mean size of L = \jL x L y = 5 nm (Fig. 2), and 10 nm (Fig. 3). The dot thickness along the z- 
direction is kept constant at L z = 3 nm, and several size ratios for each axis in the xy plane 
are shown. We use here also a value for the optical bandgap of E g = 1.51 eV. The results 
presented here were obtained using the optimized Q basis set approach discussed above. 
Notice that since this function represents all of the possible transitions of this excitonic 
system, its features would be measurable via photoluminescence excitation measurements. 
On the other hand, the photoluminescence response would correspond to the first (lower 
energy) feature in these traces, associated with the ground state of the excitonic system. 

Figure 2 shows the imaginary part of the linear optical susceptibility as a function of 
frequency for a dot with L = yj L x L y = 5 nm (a broadening of T = 2 meV is used). The 
bottom trace is for a circular dot, so that L x = L y = 5 nm. The upper two traces show 
results for elliptical dots with a size ratio L y : L x = 2 : 1 (L y = 7.0 nm, and L x = 3.5 nm), 
and 3 : 1 [L y = 8.7 nm, L x = 2.9 nm), respectively, all having the same mean size L = 5 nm. 
[Notice that although the peak heights are in arbitrary units, the ratio between different 
peaks or traces is real, reflecting the different dipole matrix elements involved.] One obvious 
difference among traces is that the first transition energy shifts to higher values by ~ 26 
and 69 meV, as the size ratio changes from 1:1 to 2:1 and 3:1, respectively. Notice that 
while increasing the size ratio, the exciton binding energies increase from 47 meV to 48 
and 49 meV for each ratio. It is then clear that the larger blue shifts are due mostly to 
the increasing confinement as the disk becomes more elliptical, and not due to the binding 
energy between electron and hole. In all the cases shown, the transition involving the ground 
state of the exciton is dominant, as the excited states appear only with smaller oscillator 
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strength. Notice further that for the larger length ratio, the spectrum is understandably 
sparser, as the levels associated with the narrow dimension are quickly pushed upwards in 
energy. Finally, since the elliptical dots have lower symmetry, accidental degeneracies are 
fewer and the transition peaks show nearly monotonically decreasing intensities (unlike the 
circular dot case). [Notice these all the susceptibility traces include transitions between 
different center-of-mass states up to the Nx = Ny = 10 levels. Additional Nx and Ny 
values would yield higher energy structure.] 

Figure 3 shows also the imaginary part of the susceptibility but for a dot size of 10 nm, 
and here with level broadening of 1 meV. The first transition energy in the circular quantum 
dot is 1.78 eV, while the values in the elliptical dots are 1.79 and 1.80 eV, for ratios 2 : 1 and 
3:1. The transition energies result here to be lower than in Fig. 2 since the confinement 
is not as strong, reducing the effective gap energy. In this set of curves, the first transition 
energy (involving the ground state of the exciton) is shifted upwards for elliptical dots, due 
to the increased confinement, although the shift is not as large as in Fig. 2, since the overall 
lengths are larger. The first excited state appears split because the two-fold degeneracy of 
the excited state is broken as the dot becomes elliptical. 

As dot size is increased, it is apparent that the geometrical effects are not as promi- 
nent, producing only a small shift of the spectrum of transitions. Incidentally, the onset of 
transitions for larger size dots (L ~ 100 nm) compares qualitatively well with experimental 
photoluminescence spectra in dots with similar disk geometry,0 where features appear in 
the energy range 1.73 - 1.74 eV, for dots with radius thought to be in the range 150-200 
nm. According to the experimental results, the additional confinements give the observed 
blue shift, compared to two-dimensional quantum-well exciton case. Similarly, a blue shift 
appears due to the elliptical shape, although for L « 100 nm they are only < 10 meV, as 
the size ratio increases to 2:1 and 3:1. 

As an example of the effects for different materials, Fig. 4 shows the imaginary part of 
the susceptibility of InAs quantum dots with lateral mean size of 12 nm and thickness of 2.8 
nm — having a dielectric constant e = 14.6, m e = 0.026mo, and heavy hole effective mass 
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rrihh = 0.41m . Here, the energy gap is taken as 0.43 eV, and we also use the optimized 
Q basis set to obtain the results shown. The convergent optimized size, Lq = 10.5 nm, is 
close to the two dimensional effective Bohr radius (a B D = 10.54 nm). In this case, the first 
transition energy shifts to higher values by « 11 and 28 meV, as the size ratio changes from 
1:1 to 2:1 and 3:1, respectively. Notice the similarity with Fig. 2, although the energy scale 
and size (12 nm) are completely different here. This similarity is due to the scaling of the 
problem in terms of a* B . For these InAs parameters, a* B = 29.8 nm, so that a* B /L xs 2.5, 
comparable to the value in Fig. 2 for GaAs, where a* B = 12.2, L = 5, and a* B /L 2.4. 

IV. CONCLUSIONS 

We have demonstrated that strong geometrical confinement effects appear on excitons in 
GaAs and InAs quantum dots with elliptical cross sections. The solutions have been obtained 
using sufficiently large basis sets as well as with an optimized basis set. The results obtained 
with the optimized basis sets, such as exciton binding energy and normalized electron-hole 
separation, are extremely close to the best converged results, and only with a relatively 
modest computational effort. 

The linear optical susceptibilities are calculated for several different lateral size ratios of 
each axis (x and y). Strong blue shifts in the susceptibilities are observed as the size ratio 
is increased and the shifts due to the geometrical shape effects are especially important for 
the smaller dot sizes (< 25 nm). The shifts are due mostly to the increasing confinement as 
the dot becomes more elliptical, and not due to the interaction energy between electron and 
hole. A splitting of the first few excited states appears in the elliptical cross section cases 
since the symmetry-related degeneracy of the excited states in the circular dot is broken. 
This gives also rise to a more monotonic decrease of the peak intensities seen as the energy 
of the transition increases. 
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APPENDIX: 

The Coulomb interaction matrix elements of the electron-hole pair can be calculated 
analytically using the harmonic-oscillator basis sets we use. With the basis set in Eq. (4), 
one can write 
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Here, the integral factor, 
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can be transformed to polar coordinates, so that, 
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where a = {iu x /h, b = fiuj y /h, I = n x > + n x — 2a — 2/3, m = n y / + n y — 27 — 25, and 
n = |(Z + m+ 1). 

The interaction matrix elements can be expressed in terms of hypergeometric functions 

as 



(n x , Tiy I 



e^/x 2 + y 



?\n x ,n y ) = ^^(T^yn'JnJn'J.nyl)- 1 / 2 



K/2] [n x /2] K/2] K/2 ] 



X £ 2 § S «!K-2a)!/5!K-2/5)! 7 !K-2 7 )!5!K-25)! 



x 



-(s^ + 1) - 7 - 5, ~(s x + s y + l)-r); ~(s x + s y ) - r] + 1; 1 - ^ 

A A A L*Jx 



x r 



-(sz + s^ + l) 



r 



^ + 1) - a - /3 



s y + 1) - 7 - 5 



/r 



(A5) 



where F is a hypergeometric function,0 r/ = a + f3 + 'y + 5, and the extra constraints of 
s y = n y + n' y = even, and s x = n x + n^. = even, are required for this matrix element to be 
non-zero. Use of this equation in the calculation of the Coulomb matrix elements was very 
important in the solution of the problem, as it reduces the computation time substantially. 
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FIGURES 

FIG. 1. GaAs quantum dot. (a) Exciton binding energy; (b) normalized electron-hole separa- 
tion as function of the quantum dot size for elliptical quantum dots with several size ratios in the 
xy plane. L x : L y shown are 1:1 to 2:1, and 3:1, with diamonds, pluses, and triangles, respectively. 
Insets: Solid, dot, dashed, and dashed-dot curves are the results for circular case and M = 30, 
100, and 500 basis sets, respectively. Diamonds show results for the M = 400 basis set with 
optimized frequency, while triangles are for M = 400 in the original u basis. Inset in (a) also shows 
(scale on the right) the value of Eq = HQ used in the calculation. 

FIG. 2. Imaginary part of linear optical susceptibility for a GaAs elliptic quantum dot with the 
same area but different axis ratios, as a function of frequency. Bottom trace is for a circular dot, 
L x = L y = 5 nm. Upper two traces show results for elliptical dots with size ratios L y : L x = 2 : 1 
(L y = 7.0 nm and L x = 3.5 nm), and 3:1 (L y = 8.7 nm and L x = 2.9 nm), respectively. 

FIG. 3. Same as Fig. 2 but for a size L = 10 nm. 

FIG. 4. Imaginary part of linear optical susceptibility for InAs elliptic quantum dots with 
lateral mean size 12 nm and thickness 2.8 nm (with T = 2.0 meV), as a function of frequency. The 
bottom trace is for a circular dot with L x = L y = 12 nm. The upper two traces show results for 
elliptical dots with a size ratio L y : L x = 2 : 1 (L y = 17 nm and L x = 8.5 nm), and 3:1 (L y = 20.7 
nm and L x = 6.9 nm), respectively. Notice the similarity with Fig. 2. 
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